Switching superconductivity in S/F bilayers by multiple-domain structures 
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We consider the effect of a multiple magnetic domain structure in a superconductor/ferromagnet 
bilayer, modeled by a ferromagnetic layer with a rotating magnetic moment. The domain walls in 
this model are of equal size as the domains, and are of Neel type. We study the superconducting 
critical temperature as a function of the rotation wavelength of the magnetic moment. The critical 
temperature of the bilayer is found to be always enhanced by the domain structure, and exhibits 
an interesting reentrant behavior. We suggest that this effect can be used for a new device where 
superconductivity may be controlled by the domain structure of the magnetic layer. 

PACS numbers: 74.45.+C, 74.62.-c, 74.78.-w 
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The study of the proximity effect between a super- 
conductor (S) and a ferromagnet (F) is currently a very 
active field due to its relevance for nanoelectronic ap- 
plications and due to the perspectives to discover new 
interesting physical phenomena. It is well-known that 
the exchange field J of the ferromagnet tends to break 
the Cooper pairs formed by electrons with opposite spins 
by acting on the electronic spins via the exchange in- 
teraction. When the exchange field exceeds the typical 
superconducting low-energy scales, J > 2ttTco (TcO is the 
superconducting critical temperature in absence of the 
proximity effect), a new, shorter length scale competes 
with the superconducting coherence length scale. For 
this case, within the quasi-classical theory of supercon- 
ductivity for diffusive structures (Usadel equations 0), 
the superconducting pair correlations have been found to 
penetrate (and also oscillate) in the F part over a short 
length = y/Df/J instead of the superconducting co- 
herence length ff = Df /2ttTco (Df is the diffusion 
constant in F) [3. Moreover, in S/F bilayers the crit- 
ical temperature Tc is suppressed by the proximity of 
the F layer and exhibits a non-monotonic dependence on 
the thickness of the F layer, df (see and references 
therein). In F/S/F sandwiches, it has been found that 
the pair-breaking effect depends on the relative orienta- 
tion of the F moments. This property led to the proposal 
of a superconducting switch [j, |5j operated by reversing 
the moment in a F layer. 

In this Letter, we introduce a new aspect of this fea- 
ture. We investigate the sensitivity of the Cooper pairs 
on the directional changes of the exchange field by study- 
ing a S/F bilayer with a moment rotating with a constant 
velocity in the F layer (see Fig. 1). This model simu- 
lates a multi-domain structure with domain walls of the 
Neel type, where the domain wall width is comparable in 
size to the size of the domains. We find that the critical 
temperature Tc is enhanced compared to that of the S /F 
bilayer with homogeneous magnetization and exhibits an 
interesting reentrant behavior. We suggest to use this 
effect for new devices where superconductivity may be 
switched on and off by controlling the magnetic domain 
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FIG. 1: Model: the moment J rotates in the F layer, 
structure. 

Recently, the theoretical understanding of the inter- 
play between the superconductivity and inhomogeneous 
magnetism has been improved by realizing that a local 
inhomogeneous exchange field induces superconducting 
triplet correlations, which penetrate over a long scale, 
in a diffusive ferromagnet 0, 0- As a result, these 
correlations lead to a significant increase of the conduc- 
tance of the ferromagnet below Tc 0. Volkov et al. 
0, 0| have demonstrated the existence of these long- 
range triplet components (TC) also in the multilayered 
F/S/F structures when the moments in each F layers are 
non-coUinear. As already emphasized in Ref. |7l. l9l llCll|. 
these TC have the property to be odd in frequency, a nec- 
essary condition for their existence in the diffusive limit. 
The study of Tc as a function of the angle between the 
moments in the F/S/F trilayer has been reported very 
recently order to clarify the situation and to 

compare with our model, we first derive generally the Us- 
adel equations in the paired-spin basis and show that the 
coexistence of singlet components (SC) with odd-triplet 
components is in fact the hallmark of any diffusive S/F 
structure 0. Short-range TC are necessarily present. 
Contrarily, the presence of TC with long-range scales re- 
quires specific conditions. For example, the long-range 
TC are shown to be absent within our model. The sen- 
sitivity of Tc on the magnetic inhomogeneity is in our 
model mainly provided by the short-range TC. 

We study the S/F bilayer within the framework of 
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the Nambu-Gor'kov formahsm and consider the diffusive 
hmit which is relevant for short mean free paths. The 
Usadel Green's function ^(R, u;„) depends on the spa- 
tial coordinate R and on the Matsubara frequency = 
TTT{2n+l) (where T is the temperature), and results from 
the momentum average of the quasi-classical Green's 
function. It is represented by a 4 x 4 Nambu-Gor'kov 
matrix in combined particle-hole and spin spaces. In 
particle-hole space it is written as 
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where the 2x2 spin-matrices g and / represent the nor- 
mal and anomalous Green's functions, g obeys a nor- 
malization condition and a nonlinear transport equation, 
coupling its components / and g Near Tc, / is small 
and g deviates slightly from its value (— i7rsgna;„) in the 
normal state, so that the Usadel transport Eq. can be 
linearized with the help of the normalization condition, 
and yields an equation for the 2x2 spin-matrix / 

DVV - »sgn(c^„) {J (R) • af + /J (R) • a*} 

-2|c^„|/ + 2^A = (2) 

with (T = [ux, (Jy, Uz) St, vector whose components are the 
Pauli matrices. In the F layer the superconducting order 
parameter A vanishes, and in the S layer the exchange 
field J is zero. The diffusion constants are D = Z^s in 
S and D = Df in F. We consider a BCS singlet super- 
conductor for which the order parameter in spin-space 
is A = As'icTy. Quite generally, / is expected to have a 
singlet and a triplet part |0| 



= fsicTy + i (ft • (t) ay 



(3) 



where ft = {ftx, fty, ftz) is the triplet vector (the basis 
formed by the vectors s = (|ti) - |iT))/2, = (|ii) - 
|tT))/2, ty = (ITT) + Ui))/2z, and = (|Ti> + |iT))/2 
is the natural basis for studying pairing states). The 
different components of / obey the system of coupled 
equations 

(DV2-2KI) = -27rA,, + 2isgn(c^„)J-ft (4) 
(DV2-2KI) ft= 2isgn(c^„)J/,. (5) 

The couplings (right-hand side) in Eq. Q-Q show 
obviously that the SC fs coexists always with at least 
one TC as early as J ^ 0. We see also generally that the 
SC is even in frequency [/s(— cj„) = fs{u>n)] as it is the 
case for the BCS superconductors, while the TC are odd 
[ ft(— a;„) = — ft(a;„)]. This is a consequence of the Pauli 
principle, which for the diffusive limit imposes the rela- 
tion /q/3(R, Wn) = — //3q(R, — Wn), whcrc a and /? arc the 
spin indices. Henceforth, we only consider positive Mat- 
subara frequencies. The Eq. Q-© are supplemented by 
the self-consistent equation relating As to fs 

- - A'/l I J I 

T 



A.ln^ = 2.r^f^- 



n>0 



(6) 



and also by the boundary conditions. The general bound- 
ary conditions at the S/F interface for the diffusive case 
have been formulated by Nazarov 0| and reduce near 
Tr to 
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where pg and pf are respectively the normal-state resis- 
tivities of the S and F metals [Eq. ||7J) follows from the 
continuity of the current at the interface] , and 

iflb dz.f)F = f)s - f)F - lb = RbA/pfis (8) 

with i?b the resistance of the S/F boundary, and A its 
area. Here z denotes the distance to the S/F interface. 
Note that near Tc Nazarov's boundary conditions are 
formall y e quivalent to the ones by Kuprianov and Lu- 
kichcv [13. At the outer surfaces of the F or S layers, 
we require that the current through the boundary has to 
vanish. It is important to note that the conditions 0-® 
do not couple the different components of /. 

If J is constant in direction, wc obtain straightfor- 
wardly that the triplet vector ft || J. This configuration 
corresponds to a triplet state with a zero spin projection 
on the quantization axis defined by J. The SC and the 
TC with zero spin projection are energetically equiva- 
lent with respect to the exchange interaction, and thus 
necessarily appear together in the ferromagnet. 

In F/S/F trilayers with an arbitrary angle between the 
moments in each F layer 0, 0| , J has not a fixed di- 
rection in the structure. If the moments are homogeneous 
in each F layer, the system Q-® can be solved easily. 
In the different F layers, the spatial dependence of the 
components of / can be written as fi{z) = fiCxp(kfz), 
leading to the algebraic equations 



Dfkj-2LUn -2iJ- 
-2i3 Dfkj-2Lu„ 



0. 



(9) 



The eigenvalue k'j is determined from the condition of 
zero determinant, which yields 

fc^± - r!„e7' ± 2<7' or kjo^n^q^ (10) 

with Qn = ^n/T^Tco- The corresponding eigenvectors 
have the form 



fs± 



fs 



1 

±J 



and 



fsQ 





fto 



(11) 

with the condition fto • J = 0. Thus, there can exist 
a priori two very different types of TC 0, • The SC 
fs and the TC || J characterized by k^^j^ have the usual 
short-range decay length in the F layer for strong J. 
On the contrary, the TC ff _L J (characterized by fc^g) 
penetrate further in the ferromagnet. Since they corre- 
spond (locally) to triplet Cooper pairs with a non-zero 
spin projection (equal-spin pairing), they are not locally 
limited by the paramagnetic interaction with J. 
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Now, we consider the S/F bilayer shown in Fig. 1 
where J rotates with a constant velocity Q in the F layer, 
i.e. J(y) = J(cos(5?/,sin(52/,0) (henceforth z designates 
no longer the quantization axis). It is straightforward to 
see that here ftz = 0. The y dependence of J is elimi- 
nated in the right-hand side of the Eq. 0)-® by con- 
sidering the new components /+ = {— ftx+ifty)^^^^ and 
/_ = {ftx + ifty)e.~^'^^ ■ The system of Eq. to solve is 

{DV^ - 2w„) = -27rA, + i J {f- - /+) (12) 
{DV^ T 2iDQdy ~ DQ^ - 2un) f±^T2i Jfs- (13) 

Since the structure is periodical in the y direction, the 
three components /; (Z = s, ±) and Ag are expanded into 
Fourier series. We note that the different Fourier compo- 
nents (labelled by p) are neither mixed by the full system 
of linearized equations nor by the boundary conditions. 
Thus, each harmonic taken separately (and determining 
a particular y dependence) represents a possible solution. 
The harmonic p realized physically is the one which gives 
the highest Tc- As a result of our calculations 0|, the 
harmonic p = yields the highest Tc. The physical rea- 
son for this is that the other harmonic solutions corre- 
spond to an inhomogencous bulk superconducting state, 
which decreases unavoidably Tc- We present here the 
equations for the case p = 0. 

In the F layer, using the boundary condition at the 
outer surface {z = —df), the components /; are thus 
sought under the form fi{z) = fi cosh [kf{z + dj)], which 
is substituted in the set of Eq. (|12|I - H13|) . This leads to 
the following linear system 

-2iC/ A-g2 /- =0 (14) 

V 2iS,f X-Q'J\f+J 

with the eigenvalue X = — QnS.j'^- This system yields 
the three eigenvalues (and consequently three fc/) 

Ae = + e^Q^-l&C/^ /2, Ao = (15) 

and the associated eigenvectors (here e = ±1) 

2!;^ )' (fco)"(j)' 

(16) 

The eigenvalues Ae give short-range decay lengths for the 
TC in strong ferromagnets, while the eigenvalue Ag de- 
termines a long-range penetration in the ferromagnet (at 
least at small Q). The solution of Eq. (|12() - (|13|1 indepen- 
dent of y is thus written as 

fi{z)= ^ aj7ijC0sh[fc/j(z-t-d/)] (17) 

where the three coefficients aj have to be determined with 
the help of the boundary conditions at the S/F interface. 
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FIG. 2: Tc versus dj for various Q(,s- Here ds = 7 nm. For 
df ~ 5 nm, corresponding to the dashed line, a switching 
effect can be observed between Q — (single domain F) and 
QCs > 0.56. 

In the S layer, the solutions for the TC f± sat- 
isfying the boundary condition at the outer surface 
(at z — ds) arc straightforwardly derived f±(jjTz) = 
c ± cosh [ks{z - ds)] where ks = y/fl„is^ + Q^, (,s = 
■\/ Ds/2'kTco, and c± are coefficients. As a result of the 
form of the eigenvector associated with the eigenvalue Aq, 
the boundary conditions for the TC at the S/F interface 
require ao = 0. This absence of the long-range TC in 
the present model is characterized by the fact that the 
triplet vector ft(y, z) is locally parallel to J(y). 

The remaining equation for the SC fs{z) can gener- 
ally not be solved analytically in the S layer since it is 
coupled with the self-consistent gap As(z) by Eq. 
We use a technical high-energy cutoff of lj„ = 21ttTcq for 
the gap equation. We iterate Eqs. l(T^ and in S nu- 
merically, together with the boundary conditions 0-©, 
which determine the coefficients Uj and c±. 

In Fig. 2 we show the dependence of Tc on df for 
different values of the dimensionless parameter Q£^s (for 
dcfinitcness we took the same parameters as in Fig. 2 
of Ref. 0: Tco = 7 K, 7 = 0.15, jb = 0.3, J = 130 K, 

= 8.9 nm, and ^/ = 7.6 nm). Tc clearly increases with 
Q. Moreover, we see that the non-monotonic behavior of 
Tc{df), responsible for the absence of superconductivity 
in a finite intermediate range of d/ (around c?/ = 5 nm in 
Fig. 2), tends to be suppressed in favor of a monotonic 
behavior. This tendency is accompanied by a reentrance 
of superconductivity as a function of magnetic inhomo- 
gcneity in the df range where Tc is zero. Fixing df in 
this range we predict that when Q exceeds some critical 
value Qcrit (~ 0.56/^^ for = 5 nm), superconductivity 
is recovered. By controlling the degree of magnetic inho- 
mogeneity it is possible to switch the bilayer between the 
normal state and the superconducting state. 

The critical behavior is illustrated for df = 5 nm in 
Fig. 3, where Tc is plotted as a function of for dif- 
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FIG. 3; Tc versus Q^s for dj = 5 nm and various ds- The 
thick curve corresponds to a cut along the dashed line in Fig. 
2. The onsets for ds < 7.44 nm behave like \/Q — Qcrit- 

ferent ds- It is clear that Tc is enhanced for any Q com- 
pared to the homogeneous case Q = 0. The switching 
behavior can be observed for superconducting layers with 
ds < dg^crit, where for the used parameters dg crit ~ 7.44 
nm. For d > ds^crit we observe for small Q a behavior 
Tc ^ Q^. For d < ds^crit there is an onset of super- 
conductivity at Qcrit, and near the onset the functional 
behavior of Tc as function of Q is Tc ~ — Qcrit- By 
choosing the thickness of the superconducting layer ap- 
propriately, it is possible to optimize the superconducting 
switch for the corresponding range of Q in which it can 
be controlled. 

Finally, we discuss the implications of our results for 
the currently debated general issue of the influence on 
the superconductivity of domain walls present in the F 
layer. The evolution of J in the domain walls can typi- 
cally be described by two kinds of models. In bulk ma- 
terial, J rotates about an axis perpendicular to the wall, 
i.e. it remains in the plane of the domain wall (Bloch 
wall). On the contrary, in thin F layers J is expected to 
evolve in the surface plane (Ncel wall). Our model can be 
seen as a model for a F layer containing many Neel walls 
whose widths are comparable to the widths of the do- 
mains. Recently, a dependence of Tc on the (controlled) 
domain state of the fcrromagnct has been observed [13 
in S/F bilayers. In the presence of many domain walls, 
Tc has been found to be slightly enhanced compared to 
the Tc obtained in the absence of domain walls. The 
authors of the Ref. ^3 argue that the domain walls in 
their experiment are likely of the Ncel type. Indeed, the 
presence of Bloch walls in a F layer induces unavoidably 
a non-zero component for J in the direction perpendic- 
ular to the layer. As a result, the exchange field acting 
on the electron spins produces in addition an electromag- 
netic field acting via the electronic charges on the orbital 
motion of electrons and thus suppresses the superconduc- 
tivity. This rising orbital pair-breaking effect competes 



with the weakened paramagnetic pair-breaking effect in 
the domain walls and is expected to reduce or even pre- 
vent any dependence of Tc on the domain state. Our 
results are thus consistent with a Tc enhancement due to 
Neel domain walls. A future step is to take into account 
more realistically the alternation of domain regions with 
domain walls. The difficulty is that the problem is intrin- 
sically two-dimensional, and the solutions with separated 
y and z dependences are no more physically acceptable 
for boundary reasons. This prevents any simple analyti- 
cal treatment even in the F layer as done here. 

The possibility to rcversibly create and remove Ncel 
walls in a controlled way as needed for the realization of 
the suggested device has been shown very recently [lij . 
The method used in Ref. is based on the exchange 
bias effect between a ferromagnetic Co thin film and an 
antiferromagnetic insulating CoO film. 

In conclusion, we have presented a model for a multi- 
ple magnetic domain structure in an S/F bilayer and have 
shown that the superconducting critical temperature is 
enhanced by the magnetic inhomogeneity and even ex- 
hibits reentrant superconducting behavior appealing for 
future applications. We suggest to use this effect for a su- 
perconducting switch operated by controlling the degree 
of inhomogeneity in the ferromagnetic layer. 
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